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Abstract - This paper explains the biophysical principles which, according to us, govern the Proteus 
mirabilis swarm phenomenon. Then, this explanation is translated into a mathematical model, essentially 
based on partial differential equations. This model is then implemented using numerical methods of the 
finite volume type in order to make simulations. The simulations show most of the characteristics which are 
observed in situ and in particular the terrace generation. 

Keywords - Proteus mirabilis swarm, modelling, partial differential equations, finite volumes. 



Contents 

1 Introduction [l| 

2 Explanation of the principles governing the swarm Q 

3 Translation of the principles into a model @ 

4 Implementation @ 

5 Simulations |^ 

5.1 Simulation 1 13 

5.2 Simulation 2 14 

5.3 Simulation 3 [19 

6 Conclusion 21 



1 Introduction 

Proteus mirabilis is a pathogenic bacterium of the urinary tract. Under specific conditions, within 
a Proteus mirabilis colony grown on a solid substratum, some bacteria, the standard form of which 
is a short cell, called "swimmer" or "vegetative" , undergo differentiation into elongated cells, called 
"swarmer" cells, capable of translocation. Having this translocation ability, those swarmer cells may 
begin to colonize the surrounding medium. After a given period of time, the swarm phenomenon 
stops and the swarmer cells de-differentiate to produce vegetative cells which would, later on, be 
able to differentiate again and resume colonization. 
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Figure 1: Proteus mirabilis: Short vegetative cells (A); swarmer cells (B); periphery of a colony (C) 
and a swarm colony with terraces (D). Scales : (A): bar — 0.5 /im; (B): bar = 5 fim; (C): bar — 10 
fj,m and (D): colony radius = 3 cm. 

This swarm phenomenon is illustrated in Figures [T] and [21 In Figure [U the Proteus mirabilis 
cells can be seen in their standard vegetative forms (A). They are about 2 /im long. The elongated 
swarmer cells which range between 10 to 25 /im and which are covered by hundreds of flagella can 
also be seen (B). The third panel (C) displays a close up of the peripheral edge of a bacterial colony 
on an agar medium filling a Petri dish. The area around the edge is the theater of fast motion of 
swarmer cell flagella making it moves while the medium is gradually colonized. The last panel (D) 
shows a colony after several hours of swarming. It forms a thin film showing the specific pattern: 
concentric rings called terraces. Those rings are made by variations in the colony's thickness. 

Figure [2] which displays three snapshots of the colony expansion on agar during an active mi- 
gration phase. The expansion process is clearly seen from the progression of the colony's edge. The 
central inoculum is located in the bottom right of each picture. 

For a precise and synthetic description of Proteus mirabilis swarm phenomenon, we refer to 
Rauprich et al. [3S] and Matsuyama et al. [33]. Among all aspects described in those references, 
periodic alternation of swarm phases and of consolidation phases, without any motion, seem to be the 
most typical one. This alternation gives rise to terraces as a resulting effect. The second one concerns 
the cell size distribution and was pointed out in Matsuyama et al. [33] . As for space distribution at 
a given time, the swarming front is essentially constituted by long cells and the proportion of shorter 
cells grows while moving towards the colony interior. Concerning time distribution at a given place, 
the colonization starts with long hyperflagellated cells and then the proportion of shorter cells grows. 
The third typical fact, pointed out in Lahaye et al. [M] [5S] is that the consolidation phase begins 
even though swarmer cells are still present in the colony front. 



2 




Figure 2: Active swarming phase on agar. The figure shows three snapshots 30 min apart during a 
migration phase on agar. 

Proteus mirabilis swarm phenomenon has fascinated micro-biologists and other scientists for over 
a century, since Hauser first described it. This fascination certainly comes from the fact that 
bacteria are often considered as very primitive organisms and from the astonishing fact that such 
primitive living forms are able to generate such nice regular structures. This phenomenon certainly 
played an important role in morphogenesis and fractal theories (see Thorn [51], Mandelbrot [30] . 
Falconer [14j and references in them) that highlighted in a very resounding way that very simple 
principles may generate very complex shapes. On the other hand, as explained by Kaufmann [22j . 
phenomena of the Proteus mirabilis swarm kind seem to belong to an intermediary step between 
unicellular organisms and multicellular organisms. As a mater of fact, key elements to the under- 
standing of evolution of life can certainly be found within the understanding of those phenomena. 

Collective behaviour of bacteria colonies has been widely observed, described and studied. For 
instance, Hoeniger [21j delineated the Proteus mirabilis swarm phenomenon and Budrene & Berg 
[To] pointed out the ability that bacteria colonies have to draw complex patterns. Shapiro [42l[43ll44] 
and more recently Rauprich et at }38| and Matsuyama et al. [33j made a wide observation campaign 
to precisely analyze the phenomenon. 

This research effort pointed out that the swarm phenomenon has to be seen as an emergent 
property leading to the consideration of bacteria as multicellular organisms (see Shapiro & Dworkin 
gni, Shapiro [331 EH] and Dworkin [H]). 

Research has mainly tackled the problems of origin and of regulation of bacterial collective 
behaviour with a genetic or a bio-chemical slant (see Shapiro [40l |4T1 l46l |47] , Belas [7] , Manos & 
Belas [3lj, Schneider et al. [39] and Gue et al. 17 \) . 

The goal of the present paper is to add a bio-physical slant concerning the swarm phenomenon. 
More precisely, we exhibit the main physical phenomena of the bacteria to achieve their colonization. 
What we have chosen to do is to make assumptions about the physical principles governing the bac- 
teria behaviour and controlling the swarm and to translate those assumptions into a mathematical 
model involving partial differential equations. Then, this model is implemented in Matlab®7 lan- 
guage. Using the resulting program, we make several simulations with several values of parameters 
involved in the model. We finally observe that macroscopic colony behaviours may be recovered 
by the simulations leading to the conclusion that our assumptions are not entirely unreasonable. 
Beyond this conclusion, the scientific approach consisting in using modelling and simulation in or- 
der to test assumptions, is another important contribution of this paper. We plan to apply it, in 
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forthcoming works, in order to refine Proteus mirabilis swarm explanation. 

Using modelling, mathematical analysis and simulation to describe or understand biological 
systems is not unusual nowadays. For an introduction to these topics, we refer to Murray 537] . Those 
kinds of mathematical methods were applied in the context of collective behaviour of bacteria. For 
a description of this, we refer to Matsuschita ^35] and Ben-Jacob & Cohen J51. 

Esipov & Shapiro |T5] set out a kinetic model describing the swarm of Proteus mirabilis at the 
colony scale recovering the main behaviours of a real colony: time alternation of swarming phases 
and consolidation phases, time alternation of colony radius increases and stops and terraces. The key 
ingredient causing their model to produce good results is the introduction of an age structuring of the 
bacterial population as is suggested in Gurtin [18] and Gurtin & Mac Gamy 19J. This introduction 
allows one to take into consideration the age of the swarmer cells when looking at their behaviour. 
This model was studied from mathematical and numerical analysis viewpoints in Esipov & Shapiro 
p] . Medvedev, Kapper & Koppel [35], Ayati [Tl[2l[l|4], Ayati & Dupont [5], Frenod [16], Ayati 
& Klapper [6l and Laurengot & Walker [26l[27|. In particular, Ayati \T showed the importance of 
using age structured model to reproduce the cell size distribution pointed out by Matsuyama et al. 

m- 

The existing models describing colonization by bacteria generally contain an advection or dif- 
fusion term allowing one to define the direction in which the colony expands conformal to reality. 
This makes the model behave in a way which conforms to reality. But the factors in front of these 
advection and diffusion terms which control the beginning, velocity, and ending of the swarm are 
mainly set from phcnomcnological considerations. 

As for the research of physical principles explaining bacterial collective behaviour we mention 
Mendelson, Salhi & Li [36 . 

Acknowledgements - We thank James Shapiro and Bruce Ayati who pointed out lacks in the 
first version of the paper and suggested ways to remedy them. We also thank Joanna Ropers for 
proofreading this paper. 

2 Explanation of the principles governing the swarm 

From a long in situ observation campaign of Proteus mirabilis colonies (see Gue et al. |17j , Keirsse 
et al. [23], Lahaye et al. \24\ I25j ) and knowledge of physical phenomena arising within complex 
fluids, we are convinced that the main physical explanations of Proteus mirabilis swarm come from 
the properties of the extra cellular matrix which is the complex fluid smothering the bacteria of the 
colony. More than ten years ago, Rauprich et al. 38J hypothesised the key role of relative osmotic 
activities at the agar colony interface. More recently. Berg [9 and Chen et al. llj further supported 
hypothesis of Lahaye et al. [M] [2S] claiming that the preservation of a critical water concentration 
allowing swarming of bacteria colonies was due to an osmotic agent. 

In short, two properties seem to play an important role: the thickness of the colony and the 
water concentration within the extra cellular matrix. 

The colony thickness plays a role in controlling the growth rate of cells in their vegetative form: 
once a given thickness is reached, the cell division ends, stopping vegetative population growth. 

The extra cellular matrix is compounded with water, polysaccharides and lipids. When the 
water concentration of this matrix is low, polysaccharides and lipids auto-organize which results 
in semi-crystallinity and visco-elastic properties. When it is high, polysaccharides and lipid auto- 
organization is lost, due to efficient solvent competition with solute-solute interactions, decreasing 
thereby the visco-elastic properties of the extra cellular matrix. 

It has been well established that the translocation ability of the colony is a consequence of the 
rotation and oscillation the hundreds of flagella of a large number of swarmer cells. On the other 
hand, observations have put into light that consolidation phases begin despite numbers of swarmer 
cells are present and moving but without any effect on colony displacement. 
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Hence, we may infer that water concentration rules the swarm phenomenon in the following way. 
When it is low, swarmer cells interacting with visco-elastic properties of the extra cellular matrix 
produce colony expansion. When the matrix water concentration is high, the visco-elastic properties 
are lost. Then swarmer cell agitation, when it occurs, does not produce a swarm. 

The differentiation process may be assumed to be a kind of dysfunction in the cell division process: 
for a small proportion of vegetative bacteria, the cell division is not brought to completion. It is also 
consistent to suppose that a cell having undergone this differentiation goes on elongating, as long as 
water concentration is low, until it reaches a maximum length depending on water concentration. 
Once this maximum length is reached, the cell oscillation and rotation begin. At this stage, the 
swarmer cells become almost metabolically dormant (see for instance Stickler [5D], Mc Coy et al. 
[33] and Falkinham and Hoffman [TS]) and then do not consume extra fuel from the extra cellular 
matrix during their migration. Then we may suppose that their life duration, in this stage, is 
finite and linked to their length. After their life duration, the metabolism-frozen swarmer cells 
dc-diffcrentiate to produce new vegetative cells. 

In a given place of the colony, the evolution of water concentration of the extra cellular matrix 
is influenced by the density of bacteria, which need water to achieve metabolism, and by the dis- 
placement of swarmer cells which bring with them a part of their surrounding medium. The third 
factor influencing matrix water concentration is the water concentration in the agar medium with 
which it exchanges water and the corresponding transfer properties at the agar - colony interface. 

3 Translation of the principles into a model 

In this section, we translate the principles explained in the previous section into a model, in the 
simplest manner possible, and even in an oversimplified manner. 

The model set out in Esipov & Shapiro [13] coupled two equations. The first of them was 
an ordinary differential equation describing the evolution of the vegetative cell biomass density. 
The swarmer cells were modelled by an age-structured density which was the solution to a partial 
differential equation describing aging, but also de-differentiation and swarming, both being age 
dependent processes. 

The model we have set out in the present paper, involves three equations. The first concerns 
vegetative cells. To be able to take into account that the differentiated cells undergo an elongation 
phase, followed by a swarming phase with no metabolism, we consider an equation modelling the 
evolution of the elongating cells and another one modelling the evolution of swarmer cells having 
stopped their metabolism, both of them involving age-structured functions. 

We begin by considering the evolution of the vegetative cells. For this, we introduce the function 
Q{t, x) which is the vegetative cell biomass density at time t e [0, T) and in a given position x e M^. 
It is supposed to be the solution to the following ordinary differential equation: 

— = — — g X(^— X(^— =— j + y p{t,a,Ka,x)e ' da. (3.1) 

The time evolution of Q is the result of cell division, described by the first term on the right-hand 
side of the equation (j3.ip . and of the de-differentiation products quantified by the second term. In 
the first term, r stands for the growth rate of the cell division process, C 1j is the proportion 
of cells undergoing differentiation, and then 1 — ^ is the proportion of cells bringing their division to 
completion. The function x has the following definition x(-^) = if < and = 1 otherwise. Then 
the factor x((£ — £)/£), where £ = £{t,x) is swarm colony thickness at t in x , takes into account 
the above evoked control of colony thickness on the cell division process in the following way: when 
£ > £, the term is set to 0. The factor xiiQ ~ Q)/Q) precludes the growth of Q by cell division 
when Q is too low, i.e. when Q < Q. 
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In the second term, p{t, a, b, x), at time t G [0, T), stands for the density of swarmer cells in position 
a: £ K^, having been freezing their metabolism at an age a £ [0, +c»), for a period of time b. Since 
the life duration of the the swarmer cells in this state is supposed to be linked to their length and 
since the length of a given swarmer cell is linked to the duration during which it elongated, we set 
that the life duration of a swarmer cell having stopped metabolism at an age a is na. The growth 
rate of the biomass of the elongating cells is the same as the one of the vegetative cells. Hence the 
biomass density in position x of swarmer cells having been stopping their metabolism at an age a, 
for a period of time b is p{t, a, b, x) e"/"^ and the integral in the second term on the right-hand side 
of (|3.ip is nothing but the biomass density of cells that de-differentiate at time t. 

The biomass density of vegetative cells at time t = is given by the following initial data: 

Q{0,x) = Qoix). (3.2) 



The description of the evolution of elongating cell population involves the following age-structured 
function ({t,a,x), which stands, at time t G [0,T), for the density of elongating cells of age a £ 
[0, -f-oo) in a; S R^. The evolution of this density is due to aging and to the transformation of 
elongating cells into swarmer cells. We assume that this transformation in a given position x takes 
place according to a transformation rate n{a, H) depending on age a of the considered bacteria and 
on water concentration H = H{t, x) of the matrix. The dependence of /i with respect to H takes 
into account that the lower H is, the larger the existence duration expectancy of elongating cells is. 

The equation modelling the evolution of ^ is: 

| + (3.3) 

Since the proportion ^ of vegetative cells that differentiate generates elongating cells with age 0, 
equation (I3.1|l and (I3.3|l arc linked by the following transfer condition, 



at,0,x) = I Qit,x) x(^) x(^). (3.4) 
The elongating cell density at the initial time is given by: 

C(0, a, x) = Co(a,2^), (3.5) 



Now we turn to the modelling of the evolution of the swarmer cells which have stopped their 
metabolism. 

We consider that the velocity of swarmer cells is the consequence of the interaction of their agitation 
with the visco-elastic properties of the matrix when water concentration is low. Then we introduce 
a function c{H) which is non-zero for the small value of H and with value for the large value of H. 
Then, as the global motion of swarmer cells, when it exists, is the consequence of cell agitation it 
is reasonable to consider that their velocity is in proportion with the colony thickness gradient Vf , 
the proportion coefhcient being c{H). In other words, we consider that the swarmer cell velocity in 
a given position x is c{H{t, x))y£{t, x). 

The definition of p{t,a,b,x) defined for t e [0,T), a G [0, +oo), b G [0, Ka) and x G was already 
given. In view of this definition, if 5 is a regular set of and fei < 62 < na are two numbers, 

/ p{t, a, b, x) dbdx, (3.6) 

Jsx[bi.b2] 

is, at time t, the density of bacteria having stopped their metabolism at an age a, for a time larger 
than bi and smaller than 62 and situated in S. The time evolution of this quantity is linked to the 



6 



aging fluxes in bi and 62 and to the swarming flux on the boundary dS of S which has the following 
expression 



{c{H)V£ p) -ndhdl, (3.7) 

9Sx[fci,62] 

where n stands for vector with norm 1, orthogonal to S and pointing outwards and where dl is the 
Lebesgue measure on dS. Hence we deduce the following equation: 

d^y^ pdhdx^ 

= / p{t,a,bi,x) dx — / p(t,a,b2,x) dx ~ / [c{H)\/ £ p) ■ n dbdl . 

Js Js JasxibiM] 

(3.8) 



dt 



Using the divergence operator V- to transform flux integral into integral over S and noticing that 
V • (c(H)V£ p) = c{H)V£ ■ Vp + pV ■ {c{H)V£), we finally obtain the partial differential equation 
p satisfies: 

If + 1^ = -c(i?) V£ • Vp - p V • {c{H)\/£). (3.9) 

Since the elongating cells that stop metabolism at an age a become swarmer cells, we have the 
following transfer condition 

p{t,a,0,x) ^ p{a,H)C, (3.10) 

and the following initial condition 

p{t,Q,b,x) = po{t,b,x), (3.11) 

to provide (|3.9p with. 

The thickness £ of the colony in a given position is in direct proportion to the biomass density 
value at this point. In other words we set 

£ ^Q + M+JV, (3.12) 

where 

r+00 

M(t,x)^ I C(t,a,x) e"/'' da, (3.13) 
Jo 

is the biomass density of elongating cells. Notice that Q -\- M. is the biomass density of metabolic 
cells; and A/", defined as 

JV{t,x)= / p{t,a,b,x) e"^^ dadb, (3.14) 

Jo Jq 

is the density of metabolically dormant cells. 

We now turn to the description of water concentration evolution. We first consider the water 
concentration of agar gel near the contact surface. Let G{t,x), defined for t G [0,T), x G M^, be 
the value of this water concentration at time t and in position x. We consider that it evolves due 
to water transfers from and to the extra cellular matrix and that those transfers are proportional to 
the water concentration difference between the two media and to a colony thickness function. We 
name 74 the proportion factor. There are also water exchanges between the agar layers close to the 
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contact surface and the deeper layers. We model those exchanges by a simple ability to relax to the 
concentration value being worth 1 with a velocity 7^. In other words, we consider that G evolves 
according to the following equation: 

dG 

— = -^,Ti£)iG-H)+Ml-G), (3.15) 

where T{£) — £ ii £ < 1 and T{£) = 1 otherwise. Equation (|3.15p is equipped with the following 
initial data 

G{0,x) =Go{x), (3.16) 

giving water concentration of agar at initial time. 

Since it is easier to model water quantity transfer as water concentration transfer, we introduce 
h(t,x), defined for t S [0,T) and x G which is the water density contained within the matrix in 
position x at time t. This density is clearly proportional to the thickness of the colony times the 
water concentration, with a proportion factor rj that is linked to the part of the swarm colony the 
matrix constitutes, i.e. h = •q£H. 

If S* is a regular subset of with boundary dS equipped with the vector field n with norm 1 
orthogonal to S and pointing outwards and with the Lebesgue measure c?Z, the quantity 



S 



hdx, (3.17) 
is the water quantity contained in the swarm colony over S. This quantity evolves according to 

/ hdxj 



dt 




I -itT{£){G-H)dx- f {c{H)V£riNH) -ndl. (3.18) 

Js JdS 



Ids 

The second term of the right-hand side of this equation quantifies the amount of water the elongating 
cells use for metabolism. This quantity is proportional, with a factor a, to the biomass density of 
elongating cells. The first term quantifies the amount of water used for metabolism by vegetative 
cells. This quantity is proportional to the biomass density of vegetative cells. The proportion factor 
is a when vegetative cells undergo cell division and is a — a' when they do not undergo cell division, 
i.e. when xii^ — £)/£)x{{Q ~ Q)/Q) = 0. Since the swarmer cells are metabolically dormant, we 
consider they do not intervene in water consumption. The third term models the water exchanges 
with agar. The last term gives the flux of water carried by the swarmer cells. In it, c{H)V£ is 
the velocity of the swarmers and then of the medium they carry with them; rjAfH is the water 
quantity carried by the swarmer cells, it is proportional to the swarmer biomass and to the water 
concentration in the matrix. 

Using the divergence operator, we deduce from ()3.18|) that h is solution to 

dh ( / (Q-Q 



dt \ V V £ y Q ' ■ 

+ itT{£){G -H)- T]c(H)V£ ■ W{MH) - r]MHV ■ {c{H)\/£), (3.19) 
endowed with the following initial data 

h{0,x) ^ ho{x), (3.20) 
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Aa -> 

— ^ — I ///„ ' — hr^ — ^ 

0-1/2 03/2 afc-1/2 afc+1/2 

ai ' a2 ' afc 

Figure 3: Elongating cell age axis discretization 



— \ ' — • — \ — ' — \ — ^ 

Xl/2 X^_i/2 Xi+i/2 Xi+i/2 

Xl X2 Xt XI 

Figure 4: Position axis discretization 



describing the water density at time t ~ Q. 
From the link between H and we deduce 

dh ( d£ dH\ 
and then, from equation (j3.19p . an equation for H 

OH Hd£ if / fQ-Q\\ ,\Q 



r] £ rj£ £ £ 



, > hn(x) 

H{0,x)^Ho{x)^ (3.23) 

!]£ 



4 Implementation 

From the model just set out, we write a Biomass Preserving Finite Volume-like numerical scheme 
(see LeVeque [211 [29]) and then we implement it. 

In this section, to simplify matters, we consider the swarm of a Proteus mirabilis population 
whose distribution remains independent of the second position variable and an agar medium which 
is also independent of the second position variable. This means that the position space is considered 
nothing but a segment, say [OjXmax], of M. We introduce a time step At, and age steps Aa and A6 
with the following link Ab — Aa — At. We also define a position step Ax such that, I Ax = a;,nax, 
for an integer /. We then define, for n G N, t„ = nAt, for fc G N*, = (fc — 1/2) Aa (see Figure [3|), 
for p e N*, bp^ {p~ l/2)Ab, and, for i G {1, . . . , /}, = (i - l/2)Aa; (see FigureE]). 

For every t„, the function Q{tn, ■) is approximated by a function with a constant value, defined 
as Qf, on every interval (xi_i/2, 2^1+1/2)1 where Xi_i/2 = Xi — Ax/2 — (i — l)Ax and Xi^i/2 ~ 
Xi + Ax/2 = iAx for i G {1, . . . , /}. Consistently, the approximation of initial data (13. 2p is: 

Q? = — / Qoix)dx. (4.1) 
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)J, for 
values of 
H near 1 



|J, for 
0<H<1 



[I for 
values of 
H near 



Aiv A(H) Ad a 

Figure 5: Shapes of /x for matrix with high, medium or low water concentration. 



Equation (|3.ip is approximated by 

= Q7 + At (i^Qr x(^^) + V?) , (4.2) 



S ^ V Q 

where and X*" are approximations of the mean values of £{tn, ■) and of /q^°° p{tn, a, na, .) e"/"^ da 
over (xi_i/2, Xi_|^i/2), supposed to be available and defined in the sequel. 

The function C,{tm ■, ■) is then approached by a function C(tn, ., .), which is constant, with worth 
Ck,i> on any rectangle TZk,i = (0^-1/2,0^+1/2) x (a;i_i/2, Xi+1/2), where afc_i/2 = Ofc - Aa/2 = 
(fc — l)Aa and 0^+1/2 = flfc + Aa/2 = fcAa, such that 



C(i„, a, a;) e^^^dadx = Aa; (e''^'^'^/^ - e('=-i)'^''/^)r 



is supposed to be close to / ({tn,a, x) e^^"^ dadx, (4.3) 



ink,, 

for fc e N* and i e {1, . . . , /}. Then, the approximation of initial data p.Sp reads 



CL = .^...A./. ^^.-nA./.^- / Co(a, x) e''/^ dad:z:. (4.4) 



Ax (g^^fj/''" e{^^i)^^/''')'7 



At each time step ^ for i e {1, . . . , /} are used for the biomass transfer from vegetative cells. 

Hence we translate condition (13.41) into: 



We introduce a non increasing function A{H), ranging in the interval [A^, Ad] with A{Q) = Ad and 
A{1) = Au,, and before approximating (|3.3[) we assume that for any value of H , function /i(., H) is 
when a < A{H) and increases very quickly after this value, as illustrated in Figure [5l In other 
words, we assume that at a given value H of the water concentration, the elongating cells remain in 
there state until an age A{H) and that they become swarmer cells very soon after this age. With 
this assumption and since Aa = At, equation p.3p may be approximated by: 

Q+i = for every ie{l,...,I} and every keN* such that (fc - 1/2) Aa < A(i7f ), (4.6) 

(^'+1 = 0, for every ie {!,..., 1} and every keN* such that (fc - 1/2) Aa > yl(i?f ), (4.7) 
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where iJ" is an approximation of H near Xi at time t„ to be defined hereafter. 

We now turn to the swarmer cehs. We approximate their density p{tm ■) by a function p which 
is constant on any parallelepiped Vk.pa = (afe~i/2, afc+1/2) x (&p-i/2, ap+1/2) x (3^1-1/2, 2:^4.1/2). The 
value pi is defined as being such that Ah p"^^ ^ is the value of p on this parallelepiped. Consistent 
with this definition, initial data p. lip is approximated by 

Pk r, i = ^ , , . , ^ , , — / po(a, b, x) e"^^ dadbdx. (4.8) 

Transfer condition (|3.10p is approached by: 

PkXi= E ^M' (4.9) 

(k-l/2)Aa>A{H^) 

for every i G {1, . . . , /}, this last relation balancing (|4.7p . 

To be able to obtain a biomass preserving scheme, we chose to discretize the form p.Sp of the 
equation giving the evolution of p. We make a splitting that consists in sequentially processing 
aging, the action of which is quantified by the two first terms on the right-hand side of equation 
p.Sp . and, translocation modelled by the last term. 

The first step of the splitting that takes into account aging, consists simply in introducing a sequence 
(Pfe.p.V^) for fc e N*, p e N* and z e {1, . . . , /} defined by 

n+1/2 n /. -1 „N 

Pk,pl ==P/c,p-i,»- (4-10) 

For the second step of the splitting we need an approximation of the flux velocity c{H)\/£ at all 
points ^ {i ^ l)Aa:: which are defined for i G {1, ...,/+ 1}, and which are the interfaces 

between the discretization interval {xi_i/2t ^i+1/2) of position space. Calling V^^^^^ the value of 
this approximation, we set i £ {2, ...,/} 

^"-V2-(zj^)^(^r-^.-r), (4.11) 

that allows one to minimize retrograde polluting propagation of the flux in regions where its velocity 
is zero. We set Vj^j = and ^"^1/2 = ^7"-i/2 which makes the swarmer cells go out of the domain 
boundary once they have reached it. 

From those flux velocities, we define the swarmer fluxes FJ^pi_ii2 i^i interfaces Xi_i/2- In an 
interface where the velocity is non negative, the value of the swarmer density just before the interface 
is taken into account. In an interface where the velocity is non positive, the value of the swarmer 
density just after the interface is taken into account. In other words we set: 

T-in T/n n+1/2 .f. y^n ^ ri 

^k,pA-l/2 - ^i-l/2Pk,p,t " 1^1-1/2 S U, 

^k,pA-l/2 — ^i-1/2 Pk,p,t-1 " ^i-1/2 > ^■ 

Once the fluxes are computed, not forgetting that p is defined for values of h smaller than na and 
that swarmer cells at an age b — na de-differentiate to produce vegetative cells, the values of the 
swarmer density approximation for the next time step is given, for every i € {1, . . . , /}, by 

Pk,p,i ~Pk,p,i ~^ ^^\^k,p,i-l/2 ^k,pA+l/2)y 

for every p e N* and fc e N* such that (p - 1/2) Afe < k(/c ~ 1/2) Aa, (4.13) 
p^+\ =0, for every peN* andke N* such that {p - 1/2) Ab > n{k - 1/2) Aa, (4.14) 
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and the approximation T)^'^^ of p{tn+i, a, na, .) e"/'^ da available for the next time step, by 

(p-l/2)Afc>K(fe-l/2)Aa ^ ^ 

(4.15) 

In every formula above the approximation (f ") of the thickness £ is given, for i G {1, . . . , /}, by 

S^^Q'l + MI+Mr, (4.16) 

feSN* 

A/-"= 5^ {e''^^^- - e^'-'^^^^^y pl^^,, (4.18) 

(p-l/2)Afc<K(fc-l/2)Aa 

in accordance with ([3A2|) . ([STTS]) and ([3A4| . 

Function G'(i„, •) is approximated by a function which is constant on each interval (a;j„i/2, 2:24.1/2)- 
Equation (|3.15|) and initial data (|3.16p are approached by 



^n+i ^Qn_ At^,r(£f )(G'/ ~ i/D + Ai7rf(l - GD, (4.19) 
Ax 

Finally, we approximate the evolution of the water quantity by 



1 r^i+1/2 

G° = — I Goix) dx. (4.20) 



+ Atj,T{£r){G7 - H-) + ^ (Ff_,/2 - (4.21) 

where, ^^"1/2 being defined by (|4.1ip . 

^"i/2=^^"i/2A/^"^r if ^"l/2<0, 

^-1/2 -^^"1/2^-1^1-1 if ^"l/2>0, 

and where /i" is the value of the constant by interval approximation of h{tn, •) on {xi_i/2T^i+i/2)- 
We then get i/f +^ by 

n+l 



1 

= -^—r. (4.23) 



The scheme just built was implemented in Matlab®7 language. 



5 Simulations 

In this section we present several simulations for various values of the given parameters Those sim- 
ulations lead to the conclusion that characteristics of the Proteus mirabilis colony can be recovered 
by the model just built. 
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Figure 6: Simulation 1 - Generation of two terraces with different time steps. 
5.1 Simulation 1 

The goal of this first simulation is to show the ability of the model to generate terraces and to study 
the stability of terrace generation with respect to time step. We set the following values of the 
parameters 

C-0.1, T = l, f = l, Q = 0.2, 7t = 0.5, 7d = 0.9, 
77 = 0.5, = 1, v4<i = 3.5 K = 2, a = 0.3, a' = 0.28, 

we take as coefficient c{H) involved in swarmer cell velocity 

c{H) = 0.2 if iJ < 0.5 and if > 0.5, (5.2) 

as position domain [0,Xmax] = [0, 1.5], and as initial functions 

Qo{x) = 0.1 if X < 0.6 and otherwise, (5-3) 
Co(a, = for all a e [0, +00) and x e [0, 1.5], (5.4) 
po{a, 6, x) = for aU a e [0, +00), b e [0, +00) and x g [0, 1.5]. (5.5) 

The initial water concentration in the extra cellular matrix is set to 0.7 and the initial water con- 
centration in agar is set to 1. The position step is set to Ax = 0.15 leading to a value of / = 11 
and simulations are done with values of At set to 0.04, 0.02 and 0.01, with Aa = Ab ~ At. The 
results associated with those values of At are given in figurelH] (line 1 gives the result for A< = 0.04, 
line 2 gives the result for At = 0.02 and line 3 gives the result for At = 0.01). In each frame of this 
figure, the horizontal axis is the position axis and the vertical one gives the biomass densities. The 
solid line is the vegetative cell density Q, the dotted line represents the sum of the vegetative cell 
density and of the elongating cell biomass density Q + if the latter is non zero. The dashed line 
is the sum of the vegetative cell density, of the elongating cell biomass density and of the swarmer 
cell biomass density Q + M + J\f, if the latter is non zero. 

We see that for the three values of the time step, the dynamics of the colony and its final shape 
are similar. The more significant differences (see last column) are located near x = 1.5. Those 
differences are due to the influence of the boundary condition and certainly also to a numerical drift 
linked to the fact that this part of the colony is generated after a long term evolution. 
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Moreover, from a quantitative point of view, if £^{At = 0.01) denotes the thickness of the colony 
at time i„ and in position Xi for the simulation made with At = 0.01, we have 



{^Fi^t = 0.02) - £f{At = OM)Y = 0.0284, (5.6) 

1=1 



i^Fi^t = 0.01) - £f (At = 0.02))^ = 0.0181, (5.7) 

1=1 

suggesting that when values of time step has an order of magnitude of 0.01 or 0.02 we are not far 
from the stability of the swarm colony dynamics with respect to time step. 

5.2 Simulation 2 

The second simulation we present, consists in setting the following values of the parameters 

^ = 0.007, T-1, £^1, Q = 0.05, 7t = 0.03, 7<i = 0.07, 
r] = 0.3, = 1, Ad^ 6.3 k = 2.5, a = 0.02, a' = 0.0194, 

in taking as coefficient c{H) involved in swarmer cell velocity the function defined by (15. 2|) , in taking 



as position domain [0,a;niax] = [0,4.5], and in taking as initial functions 

Qoix) = 0.7 if X < 0.6 and otherwise, (5.9) 

Co(a,a;) = for all a e [0, +oo) and x £ [0,4.5], (5.10) 

po{a,b,x) = for aU a £ [0,+oo), b e [0, +oo) and x e [0,4.5]. (5-11) 



The initial water concentration in the extra cellular matrix is set to and the initial water concen- 
tration in agar is set to 1. 

Concerning the discretization parameters. At = Aa = Ah — 0.05, Ax — 0.15 leading to a value 
of / = 30. The simulation is provided over 3000 time steps or, in other words, over the time interval 
[0,150]. 

The beginning of the swarm phenomenon is summarized in Figure [S] This figure is divided into 
six pictures which are in turn divided into three frames. The top frame shows the biomass densities 
at the time given just above it. The horizontal axis is the position axis and the vertical one gives 
the biomass densities. The solid line is the vegetative cell density Q, the dotted line represents the 
sum of the vegetative cell density and of the elongating cell biomass density Q + A4 ii the latter is 
non zero. The dashed line is the sum of the vegetative cell density, of the elongating cell biomass 
density and of the swarmer cell biomass density Q + Ai + J\f, if the latter is non zero. As a matter 
of fact, for instance, in the second picture the amount between the solid line and the dashed line 
has to be interpreted as the elongating cell biomass density; and, in the third picture the amount 
between the solid line and the dotted line has to be interpreted as the swarmer cell biomass density. 
The middle frame gives the water concentration H, ranging between and 1, in the extra cellular 
matrix. The bottom frame shows the water concentration G in agar. 

The first picture shows the situation at time 0.45. From the initial time on, the vegetative cell 
density has increased which has induced a water consumption proportional to the vegetative cell 
density with the proportion factor being a. As a consequence, at that time, the extra cellular matrix 
is very dry. Hence, the elongating cells which are created have a long existence duration expectancy. 
At the given time, the thickness just reaches 1. Hence the cell division stops and the only way for the 
colony thickness to grow is to increase the elongating cell biomass density. This is made possible by 
the elongation of each elongating cell, the mass of which then grows, and by the fact that elongating 
cells have a large life expectancy. 



Ax 



Aa; 
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ure 7: Simulation 2 - Beginning of the swarm phenomenon: generation of the first terrace. 
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The second picture shows the situation at time 4.45. Since time 0.45, the elongating cell biomass 
density, and the colony thickness have been increasing. The water concentration in agar has de- 
creased little due to exchange with the extra cellular matrix. In the extra cellular matrix, the water 
consumed by vegetative cells is proportional to the vegetative cell density with the proportion factor 
being a — a' . The water consumed by elongating cells is proportional to the elongating cell biomass 
density with the proportion factor being a. As a result of the water consumption and exchanges, 
at this particular time, the extra cellular matrix is little hydrated. Soon after time 4.45, elongating 
cells become swarmer cells and, as the extra cellular matrix water concentration is low, the first 
swarm step begins. 

The situation at time 5.45 (see third picture) is the following: the first swarm step is taking 
place and then swarmer cells move towards the positive a;— axis. At the same time, the extra cellular 
matrix water concentration near x = 0.8, i.e. near the colony edge, is increasing faster than near the 
colony center. As a consequence, soon after time 5.45 the extra cellular matrix water concentration 
near the colony edge goes beyond the threshold H = 0.5 stopping swarmer cells, and then jamming 
the first swarm step. 

Between time 5.45 and time 13.45 (see fourth picture), the swarm colony shape has not evolved 
a great deal. This period of time has to be interpreted as the first consolidation phase. The space 
distribution of the cell size, in accordance with Matsuyama et al. [33] , shows a large proportion of 
long cells in the colony front and a smaller proportion in the colony interior. On the other hand, the 
extra cellular matrix water concentration has grown. In particular the threshold if = 0.5 has been 
overtaken in the colony center, which will prevent the swarm from going towards the colony center 
while the next swarm step begins. 

At time 17.45 the first consolidation phase ends (see fifth picture). Since time 13.45, every 
swarmer cell has been de-differentiating to produce new vegetative cells, making the proportion of 
short cells to increase gradually, always in accordance with Matsuyama et al. [33]. As the total 
biomass density (or the thickness), near the colony edge, is low the cell division process may happen 
giving rise to an increase of vegetative cell density. At the same time, as the water consumption is 
important (proportional to the vegetative cell density with proportion factor a), the water concen- 
tration of the extra cellular matrix decreases fast. 

At time 22.95 (see sixth picture), near x = 1, the extra cellular matrix water concentration is 
low. Hence the elongating cell existence duration expectancy is long. This leads to an increase in 
thickness, despite that the thickness is larger than 1. Elongating cells turn into swarmer cells which 
cannot go towards the colony center because of the high extra cellular matrix water concentration 
in this region. Hence they go towards positive a;— axis, beginning the second swarm step. 

Figure [8] depicts the generation of a terrace in the middle of the swarm phenomenon. 

The situation in the first picture (time 60.65) is very similar to the one in the sixth picture of 
Figure [7J thickness is larger than 1, thickness increases due to elongating cells, elongating cells turn 
into swarmer cells, relative dryness of the extra cellular matrix near the colony edge and relative 
wetness of the extra cellular matrix while going towards the colony center which prevent swarming 
in that direction. 

The second picture (time 63.45) shows the situation at the end of the swarm step, i.e. at the 
beginning of the consolidation phase. The characteristics of this moment are: wetness of the extra 
cellular matrix near the colony edge and in the colony center, relative dryness with increasing water 
concentration of the extra cellular matrix between the colony edge and the the colony center. 

The colony shape has not evolved much between time 63.45 and time 73.45 (third picture) but 
the extra cellular matrix water concentration has increased. At time 73.45 the de-differentiation of 
swarmer cells is just beginning and producing vegetative cells, the population of which is growing 
by cell division. Here again, in accordance with Matsuyama et al. [33], a large proportion of long 
cells is present in the colony front and gradually this proportion decreases with time. 

The water consumption is linked to the increase yield of a relatively dry extra cellular matrix 
near the colony edge, preparing the next swarm step beginning at time 75.45 (see forth picture). 
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Figure 8: Simulation 2 - Generation of a terrace in the middle of the swarm phenomenon. 
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Figure 9: Simulation 2 - End of the colonization. 




Figure 10: Simulation 2 - Summarize of the swarm phenomenon. 



In the third and forth pictures we can see that a remaining hump around x — 1.7 has been 
created during the time period covered by the pictures of Figure [8l This hump constitutes a terrace. 

This colonization goes on until the whole position domain is conquered. Figure [5] shows the 
situation at time 100.45 and close to the end of the colonization process. The relative regularity of 
the terraces can be observed in these two pictures. We can also notice that the water concentration 
in agar and in the matrix slowly goes to 1. 

Figure [TUl summarizes the swarm phenomenon over the time period [0,150] i.e. from to the 
3000*'' time step. The horizontal axis going backward shows the time (in time step units) and the 
horizontal axis going to the right shows the position. The vertical axis represents the total biomass 
or thickness Q + Ai + M . This figure illustrates the relative regularity of the phenomenon, in time 
and position. 

Moreover, it also allows one to visualize how the terraces are formed: a large hump is generated. 
Then it drains towards the colony edge but only partially, giving rise to a remaining hump which 
makes a terrace. 

Finally, the first frame of FigurelTT] shows the same object as in Figure fTOl but restrained to time 
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Figure 11: Simulation 2 - Visualization of the relative positions of terraces, matrix water concentra- 
tion maxima and agar water concentration minima. 

interval [1000 Ai, 2000 At] — [50, 100]. The second frame shows the water concentration in the extra 
cellular matrix and the third one shows the water concentration in agar. 

This figure shows that, during the terrace formation, in front of the hump, the water concentration 
in the extra cellular matrix is high; and that, at the back of the hump, the water concentration in 
agar is low. 

5.3 Simulation 3 

The third simulation we present is carried out with the same position domain as before, with: 

$ = 0,008, T = l, £^l, Q = 0.05, 7t = 0.37, 7rf = 0.13, 
7? = 0.3, A^o = 1.2, Ad = 5.5, K = 2.5, a = 0.42, a' = 0.41, 

Qo{x) = 0.2 if X < 0.6 and otherwise, (5.13) 

function Co given by (|5.10|) . po by (|5.1ip and c by (15. 2p . The initial water concentration in the extra 
cellular matrix is set to 0.8 and the initial water concentration in agar is set to 1. 

The steps are At — Aa = A6 = 0.05, Ax = 0.15 leading to a value of / = 30. The simulation is 
provided over 4400 time steps or over the time interval [0, 220]. 

Some parameters are quite different from the previous simulation. But the main significant 
difference is that here 74 is larger than 7^. This has the following consequence, as is visible in Figure 
[13 which shows a terrace generation: the water which is contained in agar is used more. Apart from 
this difference, the way a terrace is created is very similar to simulation 1. 

The first picture shows the situation at a moment when swarmer cells arc de-differentiating and 
producing vegetative cells which are growing by cellular division. This is drying the extra cellular 
matrix. 

The second picture shows a colony thickness, near the colony edge, larger than 1. Then thickness 
increases thanks to elongating cells which are turning into swarmer cells. The extra cellular matrix 
near the colony edge is dry, while towards the colony center it is wet. 

The third picture shows the situation at the beginning of the consolidation phase. At this 
moment, the extra cellular matrix near the colony edge and in the colony center is wet. Between the 
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Figure 12: Simulation 3 - Generation of a terrace. 
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Figure 13: Simulation 3 - End of the colonization. 




Figure 14: Simulation 3 - Summarize of the swarm phenomenon. 



colony edge and the the colony center the water concentration of the extra ceUular matrix increases. 

The swarm colony shape does not evolve much until the moment of the third picture but the 
extra cellular matrix water concentration increases. 

The fifth picture shows the moment when the de-differentiation begins which will soon lead to 
the situation in the sixth picture. 

The situation in the sixth picture is similar to the one in the first picture but with a colony edge 
which has moved towards the right. 

Figure [12] shows the situation at the end of the colonization. Here again, as in Figures [U] and 
[13] we may observe the relative regularity of the terraces. 



6 Conclusion 

In this paper, we add a bio-physical slant in the Proteus mirabilis swarm description. The scientific 
method we use for this is as follows. We first express the bio-physical principles on which, according 
to us, the swarm phenomenon could be based. Then, we translate those principles into a model 
which is implemented. Simulations show that the model behaviour may recover many of the aspects 
observed in situ and listed in the introduction. This leads to the conclusion that the previously 
expressed bio-physical principles are certainly the right ones to contribute to a swarm explanation. 
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Figure 15: Simulation 3 - Visualization of the relative positions of terraces, matrix water concentra- 
tion maxima and agar water concentration minima. 

Or, more modestly, those bio- physical principles are first way-marks on the path leading to a bio- 
physical understanding of the swarm phenomenon. 

However, many questions need to be tackled. From the modelling point of view, we notice that 
the model built here is dimensionless. This means that variables have no physical unit and that 
fields are of magnitude order 1. Therefore, even if the model built here can qualitatively recover the 
swarm phenomenon, we are far from a model which can give quantitative results. In order to build 
a model for quantitative simulations, we need to take into account the physical magnitudes of the 
phenomenon and to analyze physical constant ratios. 

From the mathematical point of view, it would be nice to set out an existence result for the 
solution of the model built here. Moreover, we may notice that this model is instable. In particular, 
it has been difficult to find values of the parameters that have led to simulations showing terraces. 
The stability properties, with respect to terrace generation or to short-term evolution or to long-term 
behaviour of the model is certainly a very interesting topic which may bring out biological answers. 

From the numerical point of view, the organic assumption concerning age and time steps: A6 = 
Aa = Ai, which is needed to make the numerical scheme to be biomass preserving, will certainly 
result in bad computational results if the present method is applied for operational simulations. 
Indeed, the simulation of a real Proteus mirabilis swarm experiment requires a time step which is 
much smaller than the ageing characteristic time to resolve the dynamics in space. Hence, it is 
important to find ways to reduce the induced cost. 

The first way that can be considered, consists in noticing that the biomass preservation property 
is always realised if the ageing process is done only every v time steps, for an integer j/, setting 
Ab = Aa = i/At and modifying the method in the following way. Initializing formula (|4.ip . (|4.4p . 
KE^ and are left unchanged. 

Then for values of n G i^N* , formula (|4.2p is applied removing term 2?" and replacing At by vAt 
and formula (|45|) is applied replacing At by i^At. Formula ((46)) to (irio]) and (|4?2T|) to ((433)) are 
applied in state. 

For values of n ^ i/N, formula ()4.2p is replaced by 

Q7^^ = Q? + 'D? if n e (j^N* + 1) and Q^+^ = otherwise , (6.1) 
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formula (|4.5p is not applied, formula (|4.6p and (|4.7D are replaced by 



Cfct^ — Cfc'i for every i and every k, (6.2) 



transfer condition (|4.9p is not applied, (|4.10p is replaced by 



n+l/2 n Q^ 



and (|4.13p and (|4.14p are replaced by 



pZ]: = Pfe,;f + ^iKp,-i/2 ~ Kp,+i/2),i0T every p and fc, (6.4) 

with fluxes given by (|4.12p and (|4.1ip . The definitions (|4.16p to (|4.18p remain the same (but we can 
notice that do not evolve). Finally (g^, and ([i?^ are appHed in state. 

The second way that needs to be investigated, consists in following the method set out in Ayati 
[T] and which consists in deducing from equations p.3p and p.9p evolution equations for functions 
^(i, a + t,x) and p{t,a + t,b + t, x). From those equations we will be able to build age discrete 
equations involving age discretization possibly time dependant and not requiring an organic link 
between time and age steps. 

From the biological point of view, it would be of great interest to enrich the model in order to 
take into account genetic or chemical events. In particular, it would be nice to understand how gene 
expressions may regularize terrace generation. 
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